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ABSTRACT 

We present two epochs of coordinated X-ray-optical timing observations of the black 
hole candidate Swift J1753. 5-0127 during its 2005 outburst. The first epoch in July occurred 
at outburst peak. Two consecutive nights of observations using the McDonald Observatory 
Argos camera with the Rossi X-ray Timing Explorer show a consistent correlation with an 
immediate response and an extended tail lasting ~ 5 s. The properties of the variability and 
the correlation are consistent with thermal reprocessing in an accretion disk. The shortness 
of the lag suggests a short orbital period consistent with that recently claimed. The second 
epoch in August used the VLT FORS2 HIT mode again in conjunction with RXTE. Again 
a repeatable correlation is seen between two independent subsets of the data. In this case, 
though, the cross-correlation function has an unusual structure comprising a dip followed by a 
doubl e-peak. We suggest that this may be equivalent to the dip plus single peak structure seen 
by iKanbach et alj d2001l) in XTE Jl 1 18+480 and attributed there to syn chrotron emission; a 
similar structure was seen during later activity of Swift J1753. 5-0127 by Durant et al. (2008). 
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1 INTRODUCTION 

The X-ray source Swift J1753. 5-0127 was the first transient black 
hol e candidate discovered by the Swift m ission on 2005 June 
30 jPalmer et afl|2005l : iBurrows et~aT]|20q5h. It was also detected 



at high energies by RX TE dMorgan et all l2005h and INTEGRAL 
dCadolle Bel et aljf2005h . The X-ray s ource was asso ciated with 
an optical counterpart with R — 15.8 (Halpern 2005) and a b lue 
spectrum showing Ha and Hell emission dTorres et alj|2005al lbl). 
A variable but un resolved radio counterpart was also identified by 
lFenderetal]j2005l) . 

From the outset of the outburst, X-ray and 7-ray observa- 
tions indicated that the spectrum and variability properties of 
Swift J17 53. 5-0127 were typ i cal of a black hole sy stem in the 
hard state dMorgan et alll 20051 : ICadolle Bel et alj|2005h . The pres- 
ence of a compact radio component is also consistent with this 
dFender et al.ll2005h . While canonical black hole systems often en- 
ter a soft thermally dominated state near outburst peak, there are 
a number o f other systems that have remained in the hard state 
throughout iBrockso pp et alj2004l) . 

Relatively little is known about the binary system. We identi- 
fied significa nt optical variability and sh ort timescale correlations 
with X-rays dHvnes & Mullallvl l2005albh . Interpreting the small 
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lags seen as due to thermal reprocessing in an accretion disk around 
a black hole w e suggested a relatively short period of less than 
12hrs dHvnesl l2006). This was subsequently supported by opti- 
cal photometric observations that f ound a 3.24 hr periodicit y inter- 
preted as a superhump modulation dZurita et alj2007l l2008). If this 
interpretation is confirmed then the orbital period will be slightly 
less than this making this the shortest orbital period black hole can- 
didate yet discovered. 

In this work we present two epochs of coordinated X- 
ray and optical observations at high time-resolution during the 
main outburst in 2005, one in early July and the other in 
ea rly August. Preliminary results from this study were p r esented 
by iHvnes & Mullallvl d2005ah . iHvnes & Mullallvl d2005bh. iHvnesI 
( 2006) and IHvnesI d2007l) . We note that iDurant et al.1 d2008h re- 



ported similar, but not exactly the same behavior as we see in our 
August observation during a later phase of activity in 2007. 



2 OBSERVATIONS 

2.1 SMARTS 1.3 m Andicam data 

We performed daily monitoring, where possible, using the 
SMARTS 1.3 m telescope at Cerro Tololo, Chile, for the first three 
months of the outburst until visibility became compromised. Ob- 
servations used the Andicam dual-channel optical/IR camera. Two 
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130 s V band observations were performed together with five 50 s 
dithered H band images each night. 

All data reduction used standard procedures in IRAF. Opti- 
cal data were supplied with satisfactory pipeline reduction. IR data 
from each night were median filtered to create a sky image which 
was then subtracted from each frame before shifting and co-adding 
the individual images. In both bands aperture photometry was per- 
formed relative to a comparison star in the field. 

The optical comparison star wa s calibrated rel ative to several 
standard stars in the field of T Phe jLandold[l992h . The IR com- 
parison star was calibrated relative to 2MASS photometry. In both 
cases the uncertainty in the calibration was ~ 0.04 mag. 

Our V and H lightcurves are shown in Fig. Q] We have sup- 
pleme n ted these lightcurves with additional points fromlStill et al.l 
(2005). lTorres et all ( 120053) . and lCadolle Bel et alj J2007h . 



2.2 McDonald 2.1 m Argos data 

Fast optical photometry was obtained on 2005 July 6 and 7 for ap- 
proximately 1.5 hrs each night using the Argos CCD photometer on 
the M cDonald Observatory 2.1m telescope ( Nather & Mukadam 
|2004|) . Observations on each night were performed as an uninter- 
rupted time-series of 1 s exposures using a V filter; there was no 
dead-time between exposures. Absolute timing was synced to GPS 
1 s ticks and multiple NTP servers. 

Custom IDL software written for the instrument was used to 
subtract dark current and bias and flat-field the data before per- 
forming aperture photometry of Swift J1753. 5-0127 and a brighter 
comparison star using a 2 arcsec radius aperture. 



2.3 VLT 8 m FORS2/HIT data 

Further rapid observations w ere performed o n 2005 August 9 using 
the HIT fast imaging mode dO'Brienll2007h of the FORS2 instru- 
ment on the European Southern Observatory's UT2/VLT telescope 
on Cerro Paranal, Chile. The observation were taken with the V- 
band Bessel filter. The HIT mode uses a slit to allow only a small 
fraction of the CCD to be exposed to light. It then shifts the charge 
from this exposed region at regular intervals and stores it on the 
un-exposed regions of the CCD. Once the masked region is full (af- 
ter accumulating data for 64 s), the shutter is closed and the CCD 
read out in the usual manner. The slit-width used in the observa- 
tions was 3", which projects onto 12 (binned) pixels implying a 
time-resolution of 0.5 s. The image was very stable, with an rms 
deviation of 0.14 pixels during the entire run. The image quality 
was equally stable with a mean FWHM of 2.44 pixels (0.61") and 
an rms of 0.2 pixels. The flux was extracted using a variable 1-D 
aperture along the slit on the bias subtracted images. The width of 
the aperture was set to be 5 x FWHM, as measured from a fit to the 
comparison star profile. The sky was measured in symmetric re- 
gions around this. An identical aperture was used to extract the flux 
from the comparison star. 



2.4 RXTE/PCA data 

All high time-resolution observations were coordinated with pub- 
lic RXTEfPCA observations. For this analysis we used Standard 1 
mode lightcurves with 125 ms time-resolution and all energy chan- 
nels merged into a single bin. Each lightcurve was extracted from 
final data products using standard tools. 
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Figure 1. Lightcurves of the early outburst. The upper two panels are from 
SMARTS data, the lower from RXTE/ASM monitoring. The downward 
pointing triangle shows the Swift/UVOT measurement of Still et al. (2005). 
The upward pointing triangle is from Torres et al. (2005), and square points 
are from Cadolle Bel et al. (2007). Times of coordinated observations are 
shown. 



3 THE OUTBURST LIGHTCURVE 

In Fig. Q] we show optical, IR, and X-ray lightcurves spanning the 
first three months of outburst. The source followed a fast-rise, slow 
decay pattern with a plateau seen around 60 days after the outburst 
peak. Th is is a very comm on morphology for a transient black hole 
outburst dChen et al1ll997l) . We also indicate the times of our rapid 
optical observations. The first occurred almost exactly at the out- 
burst peak, the second about a month later on the de cline, some- 
what b efor e the plateau . The o ptical observations of IZurita et al.l 
(2008) and iDurant et al.1 d2008h were obtained much later, beyond 
the coverage here. 



4 LIGHTCURVES 

4.1 Outburst Peak 

The lightcurves from July 6 and 7, near the outburst peak, are 
shown in Fig. [2] The X-ray lightcurves show no long term vari- 
ations, but considerable rapid variability (r.m.s. ~ 12 percent on 
timescales longer than 1 s compared to a statistical uncertainty 
of 2 percent). The optical data, on the other hand, show substan- 
tial lower-frequency variability. The data appear consistent with a 

1 .2 hr period, although they are insufficient to prove that these ob- 
servations are truly periodic rather than stochastic. If real, this pe- 
riodicity does not correspond to either the proposed orbital period 
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Table 1. Characteristic frequency (i/ max ) and amplitude (r) of multi- 
Lorentzian fits to X-ray PDS. 



Date 


Component 


fmax (Hz) 


r 


Jul 6 


BLN 


0.46 ± 0.04 


0.136 ±0.004 




QPO 


0.940 ± 0.014 


0.098 ± 0.006 




QP02 


2.5 ±0.2 


0.10 ±0.02 


Jul 7 


BLN 


0.67 ± 0.06 


0.141 ±0.005 




QPO 


0.685 ± 0.005 


0.089 ± 0.003 




QP02 


1.88 ± 0.06 


0.106 ±0.006 


Aug 9 


BLN 


0.64 ±0.12 


0.20 ±0.02 




QPO 


0.251 ± 0.014 


0.045 ±0.010 




QP02 


0.87 ±0.13 


0.06 ±0.02 



Figure 3. Lightcurves from our second coordinated observation. All have 
been binned to 4 s time-resolution for clarity and an arbitrary offset in time 
has been applied (while preserving the relative timing of the X-ray and op- 
tical lightcurves.) 



~ 3.2hr dZurita et al.l [20081) or the first harmoni c of it. Spurious 
appar ent periods are not unusual in outburst (e.g. iLeibowitz et al.l 
Il99lh so we do not believe the 1.2 hr 'period' is real and disregard 
it. 



4.2 Outburst Decline 

The lightcurves from August 9 are shown in Fig. [3] The source is 
clearly considerably fainter in X-rays, but considerable variability 
is still present (r.m.s. ~ 15 percent on timescales longer than 1 s). It 
is hard to directly compare these with the July optical observations 
due to the very different characteristics of the data. Both slow and 
rapid optical variations do appear to be present, however, as in the 
previous datasets. 
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Figure 4. Power density spectra of RXTE data from the three simultaneous 
observations. The fits, each of three Lorentzians plus white noise, are de- 
scribed in the text. The July 7 PDS has been scaled down by a factor of 10 
and the August 9 one by a factor 100 for clarity. 



5 POWER DENSITY SPECTRA 

To better characterise the temporal properties at each epoch we 
have calculated power density spectra (PDS) of each RXTEfPCA 
lightcurve and show them in Fig. [4] The general form of each PDS 
is band-limited noise, resembling a broken power-law, character- 
istic of the hard state. The data from July additionally includes 
a strong low-frequency quasi-periodic oscillation close to 1 Hz. 
We fit these PDS with a multi-Lorentzian decomposition follow - 
iBelloni et all J2002l) . The band-limited noise is represented by 
a zero-centred Lorentzian (i.e. one with zero centroid frequency) 
and the QPO by a peaked one (i.e. one with non-zero centroid fre- 
quency). Even when white noise is included, we find pronounced 
residuals at a few Hz, so include a second broader QPO. T his is a 
similar model to that adopted by ICadolle Bel et al.l d2007t) except 
that we allow the third Lorentzian to be a peaked QPO rather than 
a band-limited noise component. 

With these components, good fits can be obtained for all three 
nights; we summarize the characteristic frequency and amplitude 
of each component in Table. [T] The QPOs were less pronounced 
on August 9 than they were in July but with lower amplitudes the 
same model provides a good fit, and a QPO does still appear to be 
required. We did have to fix the widths of the QPOs on August 9 to 
the average values of the July observations to adequately constrain 



the fit. The most important quantity for our later correlation anal- 
ysis is the characteristic period of the primary QPO. This is 1.06 s 
for July 6, 1 .46 s for July 7, and 4.0 s for August 9. Since the Argos 
time-resolution is only 1 s, we would not expect the QPOs to be 
resolved in the July optical data and their main effect is likely to be 
indistinguishable from an excess of white noise. The 4.0 s QPO in 
August could have a much more significant effect on our analysis. 



6 CROSS CORRELATION FUNCTIONS 
6.1 Outburst Peak 

Cross correlation functions (CCFs) for the July 6 and 7 observa- 
tions near the outburst peak are shown in Fig. [5] Although corre- 
lations were not obvious in the lightcurve they are clearly present 
at high significance, with a strong response at or close to zero lag, 
and an extended tail lasting for a few seconds. 

The obvious interpretation of the correlation is that it arises 
from thermal reprocessing either in an accretion disk or on the 
surface of the companion star. The lag of the companion star re- 
sponse modulates with orbital phase with an average of approx- 
imately the binary separation in lightseconds. With the parame- 
ters discussed in Section [7] the binary separation would be about 
5 lightseconds and so the companion peak lag would modulate in 
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Figure 2. Lightcurves from our first coordinated observations. All have been binned to 4 s time-resolution for clarity and an arbitrary offset in time has been 
applied (while preserving the relative timing of the X-ray and optical lightcurves.) 



the 0-10 s range. Since our observations span 1.5 hrs each nigh t 
and the orbital period is estimated to be 3.2 hrs dZurita et al]2008l) , 
we cover about half an orbit on each night. On both nights we ob- 
serve a peak lag close to zero, so it is unlikely that the companion is 
dominating the response, which instead appears to originate from 
the accretion disk. Additionally, black hole systems are expected 
to usually have small mass ratios, and the presence of likely su- 
perhumps ( Zurita et al.l2007l) supports this expectation in this case. 
If the mass ratio is indeed small, then the companion star will sub- 
tend only a small solid angle as seen by the X-ray source, and hence 
would be expected to reprocess much less light than the accretion 
disk. 



We should be somewhat cautious, however, as similar CCFs 
were seen in the UV in XTE Jl 118+480 (although not in the 
optical), where they may be associated with jet emission. This 
seems less likely in this case, however, as the correlations are ob- 
served at the peak of an apparently normal outburst, and the opti- 
cal spectrum is b lue and can be fitted by that of an accretion disk 
( Still et al. 2005). In addition, the optical auto-correlation function 
is found to be broader than the X-ray one as expected for thermal 
reprocessing, and unlike that seen in XTE Jl 118+480. We cannot 
conclusively rule this possibility out, however, especially as the 
CCFs seen in August (at lower luminosity) do resemble those of 
XTE Jill 8+480. 



6.2 Outburst Decline 

We have similarly calculated CCFs for August 9 using VLT data. 
We do not then have two independent observations, however we 
can compute independent CCFs by using odd and even segments of 
the optical data. For these data we computed CCFs for each optical 
segment separately and co-added them. We show both the odd and 
even CCFs, and the combined CCF in Fig. [6] 

The results are somewhat puzzling. A correlation is clearly 
present, with the structure of the feature well reproduced in both in- 
dependent CCFs. The two observatories are obviously themselves 
completely independent, so instrumental effects cannot be respon- 
sible and this feature must be intrinsic to the source. The shape of 
the CCF is strange, exhibiting a dip followed by two pronounced 
peaks. This general structure is not unpr ecedented. A s i milar dip 
followed by peak structure was seen by iKanbach et alj d200ll) in 
XTE Jl 118+480, a nd later in the activity of Swift J1753.5-0127 
Dura nt et alj d2008h found a strong dip followed by a weak peak. 

The overall dip plus peak structure, and indeed the timescales 
involved, are rather similar to those seen in XTE Jl 118+480, and 
less so to the later observations of Swift J1753. 5-0127. The double 
peak has not been seen in other sources. One possible explanation is 
that this is actually related to the weak QPO seen at this epoch. We 
deduced a 4.0 s periodicity for this, and the two main CCF peaks in 
Fig.[6]are separated by 4.0 s. We therefore interpret the splitting of 
the peak as a result of the QPO aliasing the CCF structure. 
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Figure 5. Cross correlation functions for July data from the outburst peak. 
Positive lags indicate that the optical lags behind the X-rays. No detrending 
or filtering has been done to the data. 



7 MODEL TRANSFER FUNCTIONS 

If we interpret the July correlations as arising from thermal repro- 
cessing we can investigate them furthe r by comparin g them with 
predicted disk transfer functions from lO'Brien e t al. (200 J). We 
are of course at a disadvantage in not having a full set of system 
parameters. As a starting point, we can however make plausible 
assumptions. The most reasonable is that the orbi t al peri od is in- 
deed close to 3.2 hr as suggested by IZurita et alj J2008h - Others 
are largely guesses. We assume a representative black hole mass 
of 7Mq, a mass ratio of 0.1 (reasonable for a short-period system 
believed to exhibit superhumps), an inclination of 60° (the most 
probable inclination to be drawn from a random distribution), and 
a tidally truncated disk. 

We calculate th e transfer fu nction following th e methods of 
O'Bri en et al. 1 120021) . As done bv lHvnes e t al. (2003) we then con- 
volve the transfer function with the X-ray auto-correlation function 
to obtain a predicted cross-correlation function. This is more useful 
than attempting to compare with the lightcurves themselves when 
correlated features cannot be directly seen in the lightcurves and 
the correlation only emerges when averaging many features in a 
cross-correlation function. 

We show the results in Fig. [7] Note that both model CCFs 
show an oscillatory character. This is a feature of the X-ray ACF, 
and is a consequence of the QPO in the X-ray lightcurves. Overall, 
the agreement between model and observed CCFs is good, support- 
ing the interpretation of the correlation as arising from disc repro- 
cessing. 



Figure 6. Cross-correlation and auto-correlation functions for August 9 data 
from the outburst decline. In the upper panel CCFs, the upper solid line is 
computed from the full dataset, the lower dashed and dotted lines separate 
odd and even optical observations to demonstrate the repeatability of the 
structures seen. In the lower ACFs, the solid curve is the X-ray ACF which 
has been offset by +0.1 and the dotted curve is the optical ACF. The latter 
includes a large component of statistical noise at lag 0, and may also be 
affected at small lags by the blurring of the data in the time-axis introduced 
by the data acquisition method (see Section l2~3V 



8 DISCUSSION 

We have interpreted correlations seen at the outburst peak as due to 
thermal reprocessing and those on the decline as instead associated 
with synchrotron emission as inferred for XTE Jl 118+480. These 
interpretations are obviously rather speculative and uncertain. 

Thermal reprocessing at the outburst peak is eminently plausi- 
ble and the lags are in the range expected give the proposed orbital 
period of the binary. Other X-ray binaries hav e shown corre l ations 
consi stent with reprocessing in the disk (e.g. iHvnesI Jl998l. I2006L 
2007), and the spectral energy distribution at outburst peak could 
be su ccessfully modelled as due to an accretion disk dStill et all 
2005). This can be considered the most natural interpretation. 

Associating the correlations seen on the decline with syn- 
chrotron emission is more tentative. Primarily we are interpreting 
the presence of the dip in the CCF; had this not been present, we 
would not have been led to this conclusion. In XTE Jl 1 18+480, the 
presence of optical synchrotron could be associated with a very flat 
spectral energy distribution extending into the IR. This was clearly 
not the case early in the outburst of Swift J1753. 5-0127 as the 
SED was modelled as a disk. In fact our SMARTS monitoring data 
shows no evidence for a change in the shape of the SED during the 
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Figure 7. Cross correlation functions for July data from the outburst peak 
(shown as histograms) compared with those expected for disk transfer func- 
tions (shown dashed). 

period observed, with the V — H colour remaining approximately 
constant to within the accuracy of our measurements. 

The absence of a change in the optical/IR colour does call a 
synchrotron interpretation into question. Given the evolution of the 
source between the two epochs there is no reason to expect that the 
origin of the correlated variability should be different. Nonetheless, 
the CCFs clearly are different very different, and so a change must 
have occurred. Of course, if the decline data are associated with 
synchrotron emission, and the SED has not changed, we cannot 
rule out synchrotron emission at the peak as well. This seems less 
likely, as argued above, but this possibility cannot be completely 
rejected. 



9 CONCLUSIONS 

We have presented a study of the correlated X-ray and optical 
variability in the black hole candidate Swift J1753. 5-0127 at two 
epochs, one at the peak of the outburst (2005 July 6-7) and the 
other on the decline phase (2005 August 9). Data from July shows 
a smeared and lagged correlation that can be readily interpreted 
as arising from thermal reprocessing in the accretion disk. Data 
in August show a dip and peak structure similar to that seen in 
XTE Jl 118+480. Much later in the outburst Swift J1753.5-0127 
showed a strong dip followed by a weak peak {Durant et alj |2008l) 
indicating that the CCF morphology varies with the source state. If 
it arises in the same way as in XTE Jl 1 18+480, we should interpret 
it as due to optical synchrotron emission possibly associated with a 
jet, but in the absence of spectral evidence for a jet contribution in 



the optical th is remains a s omew hat unsatisfying explanation. As 
suggested by iDurant et alj ( 120081) . it is possible that the origin is 
synchrotron or cyclotron emission from above the disk rather than 
from a jet. 
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